Domain statistics in a finite Ising chain 

S. I. DenisoJ] 

Department of Mechanics and Mathematics, Sumy State University, 
2, Rimskiy-Korsakov Street, 40007 Sumy, Ukraine 

Peter Hanggjj] 
Institut fur Physik, Universitat Augsburg, 
Universitatsstrafie 1, D-86135 Augsburg, Germany 

Abstract 

We present a comprehensive study for the statistical properties of random variables that describe 
the domain structure of a finite Ising chain with nearest-neighbor exchange interactions and free 
boundary conditions. By use of extensive combinatorics we succeed in obtaining the one-variable 
probability functions for (i) the number of domain walls, (ii) the number of up domains, and (iii) the 
number of spins in an up domain. The corresponding averages and variances of these probability 
distributions are calculated and the limiting case of an infinite chain is considered. Analyzing 
the averages and the transition time between differing chain states at low temperatures, we also 
introduce a criterion of the ferromagnetic-like behavior of a finite Ising chain. The results can be 
used to characterize magnetism in monatomic metal wires and atomic-scale memory devices. 
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I. INTRODUCTION 

The Ising model, pioneered just 80 years ago has become one of the most popular and 
useful models of statistical physics. This model system itself and its numerous generaliza- 
tions found wide application for the investigation of not only physical but also for biological, 
economical, and social systems, to name only a few. The model has also been widely used to 
characterize the cooperative behaviors in these and other systems. The salient advantages of 
the Ising model are that it is generic for systems with phase transitions, it is very convenient 
to use, and, moreover, for particular cases it can be solved exactly, i.e., its partition function 
can be calculated, at least in the thermodynamic limit, without approximations. Because 
exact solutions were found only for a certain one- and two-dimensional versions of the Ising 
model , their role for statistical physics is most important. 

The ordinary one- dimensional Ising model, which is represented by an infinite chain of 
Ising spins, i.e., spins that can either be up or down, and which do interact with each other 
via the nearest-neighbor exchange interaction, does not exhibit the ferromagnetic phase 
transition at nonzero temperatures l|. This well-known result corroborates the known 
argument of Landau and Lifshitz j^j, according to which a long-range order in infinite 
one-dimensional systems with short-range interactions is absent. The problem of long- 
range ordering, which can emerge in such systems when these conditions are violated, is 
of prominent theoretical importance. Its solution for infinite Ising chains with long-range 
interactions between s pins has been the subject of a number of remarkable studies (see, e.g., 
Refs. 0,0, E,3l9, 10, 11|). 

A priori, the statistical mechanics of a finite Ising chain with only exchange (i.e., short- 
range) interaction seems not to present an interesting topic. This is so, because it does not 
exhibit macroscopic ferromagnetic order. A detailed investigation of this model is important, 
however, by the following motivating reasons. First, the domain statistics in such finite 
chains, i.e., a probability description of forming domains, domain lengths, and domain walls 
contains most valuable information on the thermal equilibrium state. To the best of our 
knowledge, these statistics have not been studied before. The main problem is that the 
domain characteristics are not ordinary thermodynamic quantities, i.e., they are not readily 
expressed through the partition function. In short, there are no conventional methods to 
extract them. Second, a finite Ising chain represents an appropriate phenomenological model 
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for describing magnetism in monatomic metal wires deposited on substrates. Indeed, as it 



has been discovered experimentally 



12j, a Co chain on Pt substrate is characterized by the 



exchange coupling (this justifies the nearest-neighbor approximation), very large magnetic 
anisotropy (this justifies the approximation of atomic magnetic moments by Ising spins), 
and long-range ferromagnetic order (this justifies the use of finite Ising chains). We do 
emphasize that in contrast to the case with infinite Ising chains for which ferromagnetic order 
is forbidden at all nonzero temperatures finite chains can exhibit the ferromagnetic-like 
behavior (see also Sec. IV). Finally, it is very likely that the one-dimensional magnets, which 
are modelle d by finite Ising chains, will have important implication for magnetic data-storage 
technology [14 1. 

In this paper, using a variety of combinatorial approaches, we investigate thoroughly the 
domain statistics in a finite Ising chain. In Sec. II, we describe the model and introduce 
the main definitions. The joint probability functions that describe the domain structure of 
a chain are calculated in Sec. Ill by the combinatorial method. In Sec. IV, we demonstrate 
that the number of domain walls in a finite Ising chain is binomially distributed. We then 
introduce a criterion of its ferromagnetic-like behavior and consider the limiting case of an 
infinite chain. In Sec. V, we derive the probability function for the number of up domains 
and calculate its average and its variance. The probability function for the number of spins 
in an up domain (i.e., domain length) and its numerical characteristics are determined in 
Sec. VI. We summarize our novel findings in Sec. VII. Some of our technical details and 
manipulations are deferred to the Appendices. 



II. MODEL AND NOTATIONS 



We consider a finite Ising chain with free boundary conditions that contains an even 
number N of Ising spins. We assume that the spins interact through the nearest-neighbor 
ferromagnetic coupling J(> 0) and the spin variables cr, (i = 1,...,N) assume only two 
values +1 and —1, respectively, corresponding to the up and down spin orientations. For a 
given spin configuration {er^}, the chain energy is written in the form 

N-l 

E N ({vi}) = -JY, ( 21 ) 

i=l 
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According to the Gibbs distribution, the probability of this configuration is given by 

W N ({a t }) = J- e -^({-*» (2 .2) 

ZjN 

where (3 denotes the inverse temperature measured in energy units, and = 
y,{,r,} exp\—0Ew(-{cri})] is the partition function of a chain. Using, e.g., the transfer ma- 
trix method 15], Zn can be evaluated exactly, yielding the well-known result 

Z 7V = 2 iV cosh iV - 1 (/3J). (2.3) 

In order to characterize the domain distribution in a chain, we introduce the number of 
up spins, s, the number of up domains, p, the number of domain walls, i.e., the number 
of up-down and down-up spin pairs, k, and the number of spins in the first up-domain, I. 
These numbers satisfy the conditions: 0<s<N,0<p< N/2, < k < N -1, < I < N. 
These numbers are not independent because, for example, if p = then s = k = I = 0, and 
if p = N/2 then s = N/2, k = N — 1, and 1 = 1. The introduced quantities are random 
due to thermal agitation, and our main objective is to calculate the one-variable probability 
functions Pw(k), P/v(p), and Pn(1) that describe in detail domain statistics in a chain of N 
Ising spins. (Notice that some features of the probability function of magnetization have 
been studied in [16].) To this end, we also introduce the four- variable probability function 
Pn{s, p, k, I) representing the joint probability that a chain is characterized by the parameters 
s, p, k, and I. Taking into account that 

E N ({(Ti}) = E N (k) = -J(N -\) + 2Jk (2.4) 

and 

W N {{<n}) = W N (k) = -L e -« fe ), (2.5) 
this probability function can be written as 

P N {s,p,k,l) = W N {k)K N {s,p,k,l), (2.6) 

where Kn(s,p, k, I) is the number of spin configurations possessing the same set of the non- 
negative integer variables s, p, k, and I. In accordance with the basic laws of probability 
theory 17(, all the one- variable probability functions can be determined by fixing one variable 
and summing Pj^(s,p, k, I) over the admissible values of all remaining variables. 
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III. JOINT PROBABILITY FUNCTIONS 



A. Number of spin configurations 

The chain states that we describe in terms of the four variables mentioned above are, in 
general, degenerate and Kn(s,p, k, I) spin configurations correspond to each of those states. 
The states with s = p = k = l = and s = I = N, p = 1, k = are characterized by 
only one spin configuration {<7j = —1} and {(7, = +1} for all i, respectively. Therefore 
K N (0, 0, 0, 0) = K N (N, 1, 0, N) = 1. For counting K N (s,p, k,l) in other cases, when 1 < 
s < N — 1, we use combinatorial methods. Within their framework, we consider an Ising 
chain with fixed s, p, k, and I as an alternate sequence of p up boxes and k — p + 1 down 
boxes in which s up spins and N — s down spins are distributed. 

Because the first up box must contain / up spins and in each other up box must be at 
least one up spin (hence the condition s — I > p — 1 must hold), the number Mf of different 
distributions of s up spins over p up boxes equals C p s Zi_i- Here, the binomial coefficients C™ 
with integers n and m are defined as follows: C™ = n\/(n — m)\m\ if n > m > 0, C™ = 
if m > n > or if n > and m < 0, and C° n = C" = 1 for all integer n. Using these 
properties, we can represent M T in the form M T = CgZf-^ apkl (A spM = 5 s ,o5 P ,o4,o^,o, K,m is 
the Kronecker symbol), which is valid for < s < N . 

Similarly, the number of different distributions of N — s down spins over k — p + 1 
down boxes, each of which contains at least one down spin, is given by Mj = C^~F S _ ± . By 
the same reason as in the previous case, this formula is also valid for all values of s. It 
may seem at a first glance, due to the multiplication principle of combinatorics, that the 
representation K N (s,p,k,l) = M^M^ is valid. This is, however, generally not the case. To 
obtain the correct formula for Kn(s,p, k, I) we note that for p > 1 the variable k can take 
only three values: k — 2p, k — 2p — 1, and k = 2p — 2. If k = 2p (k = 2p — 2) then both, the 
first and the last domains in a chain belong to the down (up) type, and the previous formula 
is valid. However, if k = 2p — 1 then those domains belong to different types and, since the 
first domain can be either of up or down type, we find for this case K N (s,p, k, I) = 2M^M^. 
Collecting the above results, we obtain 

K N (s,p, k, I) = (1 + 8 k ^ClZ^f spkl C k N -V s _ v (3.1) 

Note that this representation of the function K N (s ) p, k, I) is valid if the values of its variables 

5 



are compatible with each other. 



B. Three-variable joint probability functions 

By using Eqs. (12. 6(1 and (j3.1|) . we next can determine all the three- variable joint proba- 
bility functions, namely, P^(s,p,k), P^{p,k,l), Pn(s,p,1), and P/v(s,/c,/). In view of our 
purpose, however, i.e., for determining the mentioned above one- variable probability func- 
tions, we need only two of them, Pn(s,p, k) and Pn{p, k, I). According to the common rule, 
to calculate the joint probability function P]\r(s,p,k) we need to fix its variables and to 
perform the summation of Pn{s,p, k,l) over the admissible values of I. Since I = (and 
p = k = 0) if s = 0, and I = N (and p = 1, k = 0) if s = N, and 1 < I < s - p + 1 if 
1 < s < N - 1, we find 

f W N (0), s = 0, s = N, 
P N (s,p,k) = \ (3.2) 
[P N (s,p,k), l<s<N-l. 

Here, we have used the conditions that P N (0, 0, 0, 0) = W N (0) and P N (N, 1, 0, N) = W N (0), 
and introduced the notation 

s— p+1 

P N (s,p,k)= P N{s,P,k,l). (3.3) 
l=i 

In order to evaluate the above sum, we evaluate first the sum Si = Y%=V~ ^a-i-v ^ 
p = 1 (i.e., s > 1) then, using the properties of binomial coefficients, we obtain S\ = 
£?=i C'}^ = CZ{ = 1, and if p > 2 (i.e., s > 2) then, using the relation Q, 



n 



m=0 

being valid when the binomial coefficients exist, we have S\ = CgZ 1 — CpZ^- Since C^l\ = 
if p > 2 and C P S Z\ = 1 if p — 1, we conclude that the formula = Cf.j^ holds for all s > 1 
(we recall that 1 < p < s). Therefore, substituting Eq. (|3.1|) into Eq. (|2.6j) . from Eq. (|3.3|) 
we obtain 

P N (s,p, k) = (1 + 4, 2p -i)^(A;)C , s p r 1 1 C^_ 1 . (3.5) 

Although this formula has been derived for 1 < s < N — 1, its right-hand side exists also 
for s = (when p = k = 0) and s = N (when p = 1 and fc = 0). Furthermore, since 
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Pjv(0,0, 0) = Pn(N, 1,0) = Wn(0), the desired joint probability function ()3.2|) is given by 
the same expression, i.e., 

P N (s,p, k) = (l + S k , 2p -i)W N (k)C p -lC k N -_ p s ^. (3.6) 

To evaluate Pn(p, k, I), we need to find the admissible values of s for fixed p, k, and I. If 
p = then s = (and p = I = 0), if p = 1 then s = I for k = 1,2 and s = I = N for = 0, 
and if 2 < p < iV/2 then p + 1 — 1 < s < iV — (A; — p+ 1) (recall that — p + 1 is the number 
of down domains in a chain). According to this observation we get 



P N (p,k,l) 



W N (0), p = 0, 

P N (l,l,k,l), p=l, (3.7) 
[P N (p,k,l), 2<p<N/2, 



where 



and 



P N (l, 1, k, I) = (1 + 5 kA )W N (k)C k N zU (3.8) 



By use of the relation 



in evaluating 







N+p-k-1 

P N (p,k,l) = Yl Ms,P,k,l). (3.9) 

s=p+l-l 



^ Cc+m,Cb-m — (3.10) 
m=0 



N+p-k-1 

W-J-lViV-a-l — u JV-!-li \ 6 - LL ) 

s=p+l-l 



from Eq. ()3.9j) we obtain 

P N (p,k,l) = (1 + S k ,2 P -i)W N (k)C^ r _} l _ 1 . (3.12) 

Comparing this formula with (|3.8|) . we check that, although Eq. (|3.12j) is derived for p > 2, 
it remains also valid for p = 1. Therefore, introducing the notation = 5 P) o^fc,o^z,o 5 the 
result in Eq. ()3.7|) can be represented in the appealing form 

P N {p, k, /) = (! + 5 k ,2 P -i)W N (k)C k N _)lf pkl . (3.13) 
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C. Two- variable probability functions 



The four-variable joint probability function Pjy(s,p, k,l) generates six different two- 
variable joint probability functions. But keeping in mind the one-variable probability func- 
tions, we calculate only two of them, namely P^{p^k) and P^(k,l). Because P N (p,k) = 
Wn{0) if k = and the parameter s varies from p to N + p — k — lifl<k<N — 1, the 
former is given by 

( W N {0), k = 0, 
PN(p,k) = \ (3.14) 
(P N {p,k), l<k<N-l, 

where 

N+p-k-l 

P N {p,k)= £ P N {s,p,k). (3.15) 

s=p 

Taking into account that, according to Eq. (|3.10J) . the relation 

N+p— fc-1 

Cs-iCjv-s-i — Cjv-i (3.16) 

s=p 

holds (we used the condition C™~ m = C™), we obtain 

P N (p, k) = (1 + 4, 2p -i)^(A;)C^_ 1 . (3.17) 

The same expression for Pn{p, k) follows also from the joint probability function (J3.13j) . 
Indeed, since for 1 < k < N — 1 the parameter / is varied from 1 to iV — k (the maximal 
value of / corresponds to the case when all the remaining p — 1 up domains and all k—p+1 
down domains consist of one spin), we have 

N-k 

P N (p,k)= P N {s,p,l). (3.18) 

Substituting Eq. 1CTT3I) into Eq. (I3~T£J) . and using the formula E^Cv-Vi = Cjv-i-Cft-n 
which results from Eq. ()3.4|) . and granting the conditions C|_ x = {k ^ 0) and C%_ x = 1, 
we again arrive at Eq. ()3.17|) . 

In order to find P/v(^, from Eq. (J3.13)) . we first notice that for fixed k the parameter p 
can take only one or two values. More precisely, if k is odd, i.e., k = 2h+l (0 < h < N/2 — 1), 
then the first and the last spins of a chain belong to the different types. In this case, p = h+1, 
1 < I < N — 2h — 1, and P/v(^, = -Pjv(^ + 1,2/t + 1,Z). On the contrary, if is even, 
i.e., = 2h, then the first and the last spins belong to the same type. In accordance with 



this, at 1 < I < N — 2h, the parameter p takes two values p = h (if a chain begins and 
ends by the down spins) and p = h + 1 (if a chain begins and ends by the up spins), and so 
P N (k, I) = P N {h, 2h, I) + P N (h + 1, 2h, I). Moreover, ifZ = 0(p = Jfc = 0)orZ = JV(p = l, 
k — 0) then P N (k,l) = Wn{0). Combining these results yields 

W N (0), I = 0, I = N, 
P N (kJ) = { (3.19) 

P N {k,l), 1<1<N-1, 
where 

[fc/2]+l 

P N (kJ)= J2 PN(P,k,l) (3.20) 

p=[(fc+l)/2] 

wherein [x] denotes the integer part of x. Finally, taking into consideration the following 
relation 

[fc/2]+l 

E (l + ^,2 P -i) = 2, (3.21) 

p=[(fc+l)/2] 

we find from Eqs. (|3.19j) . (|3.20|) . and (|3.13j) for the desired probability function the result 

P N {k, I) = (2 - 5 lfi - 6 l<N )W N (k)C k N -^tt kl (3.22) 

(A^ = Skfi^i,o)- 

IV. DISTRIBUTION OF DOMAIN WALLS 

According to Eqs. (|3.17|) and ()3.2ip . the one- variable probability function 

[fc/2]+l 

P^(A;)= E P N {p,k), (4.1) 

p=[(fc+l)/2] 

which characterizes the distribution of the number of domain walls in a finite chain, is 
written as 

P N (k) = 2W N (k)C h N _ v (4.2) 

The last formula reflects the fact that Pjsr(k) is the overall probability of all 2C^__ 1 spin 
configurations each of which possesses k domain walls and has the probability Wj\r(&) (we are 
grateful to an anonymous referee for this point). By using Eqs. fl2.3|) - (j2.5j) and introducing 
the designation r = (1 + e 2/3J ) _1 (0 < r < 1/2), it can be recast to read 

P N (k)=C k N _ 1 r k (l-r) N ~ 1 - k (4.3) 



(0 < k < N — 1). This explicitly shows that a binomial distribution for k emerges. 

The fact that the number of domain walls in an Ising chain is binomially distributed has 
a simple interpretation. To demonstrate this we first remind ourselves that the binomial 
distribution gives the probability C™q m (l — g) n ~ m of m successes in a sequence of n inde- 
pendent trials, called Bernoulli trials, each of which has only one outcome, i.e., success with 
probability q or failure with probability 1 — q. In our case, we consider a chain as a result 
of one-by-one addition of iV — 1 spins to the seed one. We treat each addition as a trial 
whose outcome is either along or opposite direction of the added spin with respect to the 
direction of the nearest spin. If the added spin has the opposite direction then the domain 
wall is formed and we call such outcome a success. Hence, a chain of iV Ising spins with 
k domain walls is equivalent, in the above sense, to a sequence of iV — 1 Bernoulli trials 
that have k successes. Due to the exchange interaction, the probability q of success equals 
e~@ J / (e _/3J + e@ J ) = r, and so the probability that a chain has exactly k domain walls is 
indeed given by Eq. (|4.3j) . 

The probability function (|4.3jl is properly normalized, i.e., P w(k) = 1, and, in 

accordance with well-known properties of the binomial distribution [17], the average (k) = 
Ek=i kP N (k) and the variance o\ = (k 2 ) - (k) 2 = J2k=i k 2 P N (k) - (k) 2 of the number of 
domain walls in a chain assume the form 

(jfe) = (N - l)r, u\ = (N- 1)(1- r)r. (4.4) 

For (3 J < 1, we obtain (jfe) — (N — 1)(1 - (3J)/2 and a\ = (N - 1)(1 - (3 2 J 2 )/A with linear 
and quadratic accuracy in (3 J, respectively. The relation lim^j^ (k) = (N — l)/2 makes 
explicit that in the high-temperature case, which is characterized by the condition r = 1/2, 
approximately one domain wall falls on two spins, implying that each domain contains, on 
average, two spins (see also Sec. VI). 

The increase of j3J leads to the decrease of r and Eq. (|4.4jl yields (k) w of « (N— l)e~ 2/3J 
for large enough values of (3 J. If 2(3 J 3> In AT then the condition (k) <C 1 holds, which 
indicates that in this case the spin configurations {<jj = 1} and {(jj = — 1} play the main 
role in determining the chain properties. Let r tr and r m be the transition time between 
these states, i.e., the average time during which a chain passes from the state {<Ji = 1} to 
the state {cr, = —1}, and the measurement time, i.e., the total time necessary to perform a 
measurement of the magnetization, respectively. Then, if 2(3 J 3> In iV and Tt r 3> r m , a chain 
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possesses a spontaneous magnetization. In other words, these conditions form a criterion of 
the ferromagnetic-like behavior of a finite Ising chain. Notice that in the thermodynamic 
limit (N — > oo) the second condition holds always, see just below, while the first condition 
holds only if T = 0. Therefore, in full agreement with 1], the long-range ferromagnetic 
order in an infinite chain occurs only at zero temperature. 

In order to estimate the dependence of r tr on N, it is necessary to go beyond the Ising 
model. To this end, we consider the Ising spins as the classical Heisenberg spins with 
large uniaxial anisotropy and use the Arrhenius-Neel law l^J [^J. According to it, the 
average time r between the spin reversals can be evaluated as r ~ tqe^ 3 , where To is the 
spin precession time, and AU(^$> is the height of the potential barrier between two 
equilibrium directions of each spin. Since a chain in a state with (k) <C 1 can be associated 
with a single enlarged spin for which the potential barrier height is given by NAU, we find 
that the transition time r tr exponentially grows with iV: r tr ~ Te {N-±)P&u _ Note also that 
because (3 J — > oo and r tr — > oo as T — > 0, there is always a temperature interval where a 
finite Ising chain exhibits the ferromagnetic-like behavior. 

We briefly discuss here also the problem of domain walls distribution in an infinite chain. 
As is well known [21], the binomial distribution has no unique asymptotic as the num- 
ber of Bernoulli trials tends to infinity. However, since in our case the parameter r does 
not depend on N, the probability function P^ik) does have it. To characterize Pjsr{k) as 
N — ► oo, we assume in Eq. ()4.3j) that k — (k) + o^z and define the probability function 
Pk{z) = liniAT^oo crfeP/v((^) + &kz) of the parameter z. Applying a local limit theorem j^L], we 
immediately find that Pk(z) has the standard Gaussian distribution: Pk(z) = {2ix)~ l l 2 e~ z2 1 2 . 



V. DISTRIBUTION OF UP DOMAINS 



To derive the one- variable probability function Pn{p) that describes the distribution of 
the number of up domains in a finite chain, we again proceed from the joint probability 
function Pn(p, k). A simple consideration shows that k = if p = 0, k = 2p — i (i = 0, 1, 2) 
if 1 < p < N/2 — 1, and k — 2p — i (i — 1, 2) if p — N/2. Hence, for fixed p the parameter 
k is varied from 2p — 2 + 25 P)0 to 2p — S Pi n/2, and Pn(p) is given by 

2 P-' 5 p,JV/2 

Pn(p)= E P*(p,k). (5.1) 

fc=2p-2+2<5 Pi0 
11 



Substituting Eq. (|3.17|) into this formula and taking into account the properties of binomial 
coefficients, this probability distribution is obtained as 



Pn(p) = E(l + S ltn )W N (2p - n)C 2 ^. (5.2) 

n=0 

Finally, using Eqs. ()2.3|) - ()2.5|) . from Eq. ()5.2j) we find the desired probability function in 
the form 

Pn(p) = \ E(l + M^-iV^l - rf' 1 - 2 ^ (5.3) 

* n=0 

(0 < p < N/2). This distribution, due to its formal closeness to the ordinary binomial 
distribution, will be termed the modified binomial distribution. 
Taking into consideration the results for the finite series Q| 

[(n xT Cf +V = (1 + ^ T ~^ ~ v/i) " (5.5) 
and the properties of the binomial coefficients, one readily finds that the quantities 

N/2 

c=e ^ p rrrr 2p ~ n (i - r)^- 1 - 2 ^" (5.6) 

p=0 

(n, m = 0, 1, 2) can be represented in the form 

1 + (-l) n+m (l - 2r) N - 1 - rn ] . (5.7) 



m 

jrn ' 



2 

With these results, it follows that the modified binomial distribution is also properly nor- 
malized, namely, 

N/2 1 2 

E^W = -E( 1 + M4° = i- (s.s) 

p=0 / n=0 

The average of the number of up domains in a chain is defined as (p) = Ylp=oP Pn{p)- 
Using the probability function ()5.3|) and the identity 

2p c^_7 = nC 2 N p I? + (N- i)c 2 N p :r\ (5.9) 

which can be verified directly, we find 

<P> = 7 EC 1 + <MK + (iv - (5-10) 

4 n=0 
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and substituting Eq. ()5.7|) into Eq. ()5.1U|) we obtain 

(P) = \ + \{N-I)r. (5.11) 

It may seem strange at a first glance that (p) = 1/2 in the low-temperature limit (r — > 0), 
but in compliance with Eq. ()5.3|) at r — > only two states of a chain, namely {<7j = — 1} 
(p = 0) and {<Ji = 1} (p = 1), have nonzero probabilities and they are equal to 1/2. Notice 
also that, according to Eqs. ()4.4j) and (|5.11jl . the general condition 2(p) = 1 + (fc) always 
holds. 

To find the variance a 2 = (p 2 ) — (p) 2 of the number of up domains in a chain, we first 
calculate the second moment (p 2 ) = Y^p=oP 2 Pn{p)- By applying the identity 

4p 2 C*ri n = n 2 C%I? + (2n + 1)(N- lJC^" 1 

+ (AT-l)(iV-2)C^- 2 (5.12) 

(note that last term equals zero at iV = 2) and the notation (|5.6|) . this quantity can be 
expressed as 

8 n=0 

+ (A r-i)(iV-2)/ 2 ]. (5.13) 

Inserting Eq. (|5.7j) into this formula and performing straightforward calculations yields 

(p 2 ) = 3 - + ^(N-l)r+ 1 -(N-l)(N-2)r 2 

+ l(l-2r) iV - 1 . (5.14) 

Therefore, using the definition of the variance a 2 , we find 

*l = I + \(N - 1)(1 - r)r + i(l - 2r) N -\ (5.15) 

The fact that (p 2 ) — »■ 1/2 (cr 2 — > 1/4) as r — > has the same interpretation as the low- 
temperature behavior of (p) given above. 

In conclusion, we note that if p = (p) + a v z as N — > oo then the parameter z again 
possesses a standard Gaussian distribution (see Appendix A). 
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VI. DISTRIBUTION OF DOMAIN LENGTHS 



To find the probability function of domain lengths, Pjv(0> we proceed from the joint 
probability function (J3.22)) . Since for 1 < I < N — 1 the number of domain walls k can vary 
from 1 to N — I and k = if I = or I = N, this probability function can be written as 



Pv(0 = _ (6.1) 




where 

N-l 

P N (l) = £iVM). (6-2) 

k=l 

In virtue of this, taking into account that Wjv(O) = (1— r) N ~ 1 /2 and using the standard series 
J2k=o^n xk = (l + ^) n which permits us to reduce Eq. ()6.2|) into the form Pjy(Z) = r(l — r) 1 ^ 1 , 
we obtain 

f (1 -r) N - l /2, I = 0, I = N, 
Pn(1) = { (6.3) 
{r(l-r) l -\ 1<1<N-1. 

It is not difficult to verify, using the well-known relation, J2k=o xk = (1 — x n+1 )/(l — x), 
that this distribution, which we term the finite geometric distribution, is normalized, i.e., 
J2iLoPn(1) = 1- Note also that in the limit of an infinite chain the domain lengths distribu- 
tion ()6.3|) is reduced to the geometric distribution, Poo(Z) = r(l — r)' _1 (/ = 1, 2, . . .), whose 
mean and variance are 1/r and (1 — r)/r 2 , respectively. 

By applying the standard series, J2k=o = [ x + ( nx — n— l)x n+1 }/ (1 — x) 2 , the average 
length of an up domain, (I) = J2i=o ^Pat(O) can represented in the form 

= 2 - iVrq-r)^- 1 - 2(1 -r)^ _ 

An alternative derivation of this result is presented in the Appendix B. According to this 
expression, we find limAr^ 00 (/) = 1/r, lim r ^o(0 — N/2, and lim r _^i/2(Z) = 2 — (N + 2)2~ N . 
The last condition shows that in the high-temperature limit the average number of spins 
that form one up domain in a long chain is approximately equal to 2. 

All other moments of the finite geometric distribution (|6.3J) are also calculated exactly. 
In particular, the variance of domain lengths, of = (I 2 ) — (I) 2 , is given by 

of = (1 - r)/r 2 - (N/2) 2 \2 + (1 - r)^ 1 }^ - r) N ^ 
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- N[l - 2r + (1 - r) N ]{l - rf^/r 

+ [r-(l-r) N ](l-r) N /r 2 . (6.5) 

With this result we immediately obtain af —>■ (1 — r)/r 2 as N — > oo, of — > N 2 /A as r — > 0, 
and a 2 -> 2 - (N 2 - 2)2' N - (N + 2) 2 2~ 2Af as r -> 1/2. 

To gain more insight into the domain statistics, we also introduce the probability function 
P^(l) that describes the domain lengths distribution in assemblies of Ising chains each of 
which contains at least one up domain of nonzero length. In other words, we assume that 
the parameter I can vary from 1 to N. In this case, in contrast with Eq. (|3.22j) . the joint 
probability function P£(k,l) that a chain contains k domain walls and the first up domain 
contains I spins is written as 

P+(k, = (2 - Si,N)W+(k)C k N zU, (6.6) 

where W£(k) = e~ pEN{h) /Z^ and Z% = Z N - e ^ J(Ar_1) is the partition function for 
such a chain. Therefore, the probability function P^(l), which is defined as P£(l) = 
Y!rk=i-6i N PnQ^i 0) assumes the form 

2 [(l-r)*- 1 ^, l = N, 

1 ^ r > [r(l-ry-\ \<l<N-l. 

One can again verify that the normalization condition Ya=\Pn(J) = 1 holds, and that 
the average length of an up domain, + = J^fLi ^jv(O) can ^ e represented as 



2 - Nr(l - r) 1 *- 1 - 2(1 - r) N 

«> + = *- rl -rA ■ (6 ' 8) 



As N ^ oo, the averages (/) + and (I) tend to the same limit, but their low- and high- 
temperature limits are different: lim r ^ o (0 + = N an d nm rwi/2(0 + = 2 — N/(2 N — 1). The 
former confirms the possibility of the existence of the ferromagnetic-like state in a finite 
Ising chain at low temperatures. We note in this context that the condition 2/3 J 3> In TV, 
which we introduced in Sec. IV, is equivalent to the condition £ 3> N, where £ is the spin- 
spin correlation length. According to [3], in the case of Ising chains with only exchange 
interaction and free boundary conditions the correlation length is given by the exact formula 
£ = — ln -1 (l — 2r). Since at low temperatures the asymptotic relations £ ~ l/2r and 
r ~ e~ 2l3J hold, the condition £ iV is actually reduced to 2/3 J 3> IniV. 
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VII. CONCLUSIONS 



We have determined the domain statistics in a finite chain of Ising spins that interact 
only through the exchange interaction. For a chain with an even number of spins and free 
boundary conditions, we have calculated, via a combinatorial approach, the joint proba- 
bility function of four random variables (namely, the number of up spins, the number of 
up domains, the number of domain walls, and the number of spins in the first up-domain) 
that thoroughly describe the domain statistics in a chain. Starting out from this result, we 
derived the probability distribution functions for the number of domain walls, number of 
up domains, and number of spins in an up domain. The first corresponds to the binomial 
distribution, the second to the modified binomial distribution, and the third to the finite 
geometric distribution. For each of them, we have calculated the corresponding thermal 
average and variance, have analyzed the cases of low and high temperatures, and, as well, 
have considered the thermodynamic limit. 

In addition, we have derived a criterion that a finite Ising chain exhibits the ferromagnetic- 
like behavior. According to it, the transition time between the fully magnetized chain states 
must exceed the measuring time, and the average number of domain walls must be much 
less than 1. These conditions hold, i.e., a finite Ising chain does display a ferromagnetic-like 
order on the measuring time scale, if the temperature is sufficiently small. 
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APPENDIX A: DISTRIBUTION OF DOMAINS IN AN INFINITE CHAIN 

To find the probability function of the parameter z, P p (z) = limjv-»oo ^pPn{{p) + cr p z), 
first we represent the binomial coefficients in Eq. (J5.3|) as 





(Al) 



r(6 + l)r(a-6 + l) 
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(r(x) is the gamma function) and use the Stirling formula |23( 

~ V2^e~ x x x - 1/2 (x -> oo). (A2) 



This yields 



and 



,2{p)+2a p z-n eR n ( 1 - T\ JVr+2o- p *+3/2-r-n 



c *vn"p*-n _ » ^_ - < (A3) 



/ 2CT„2; + 1 — r — 72\-A r (l-r)+2<7 p z+ 

V ^ — ) 



X 



+ 



3/2-r-n 

iV(l-r) ' 

2a p Z + 2 — r — n\-iVr-2<Tj,2-3/2+r+n 



(A4) 



as iV — > oo, and so 



err * 

P p (z) = lim / £(1 + M^n- (A5) 

Finally, taking into account that a^/N — > r(l — r)/4 and lni? n — > —z 2 /2 — 1 as iV — > oo, 
we indeed find that P p (z) = ^Tr)" 1 ^-* 2 ^ 



APPENDIX B: ALTERNATIVE DERIVATION OF EQUATION 

Using the joint probability function (|3.6jh we can also represent (I) in the following form: 

N/2 -I 2p-<5 PjJV/2 Ar+p_fc_i 

<*> = £- E E *JM*,p,*)- (bi) 

p=l P fc=2p-2 s=p 

Since sC^Z\ = pC p , Wtf(k) = r k (l — r) N ~ k ~ 1 /2, and according to the result for the series 

E sC^lC k N ZU = pC k N +1 , (B2) 



Eq. (jBljl can be rewritten as 



d_ r \N-i 2 

(0 = - ^ E(! + M^, (B3) 

/ n=0 



where 

AT/2 



n = E i 

p=i 1 



r JP-n crl _ n (B4) 



r 
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Upon calculating these quantities with the help of the series ()5.4j) and (|5.5p . 

(n = 0, 1, 2), and substituting the corresponding expressions into Eq. (|B3|) . we again obtain 
Eq. flOp . 
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